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A measurement of the proton spin structure function gf(x, Q 2 ) in deep-inelastic scattering is 
presented. The data were taken with the 27.6 GeV longitudinally polarised positron beam at HERA 
incident on a longitudinally polarised pure hydrogen gas target internal to the storage ring. The 
kinematic range is 0.021 < x < 0.85 and 0.8 GeV 2 < Q 2 < 20 GeV 2 . The integral J°'^ gf(x) dx 
evaluated at Q 2 of 2.5 GeV 2 is 0.122 ± 0.003(stat.)±0.010(syst.). 



Deep-inelastic lepton-nucleon scattering is well estab- 
lished as a powerful tool for the investigation of nucleon 
structure. The measured structure functions have been 
successfully interpreted in terms of parton distributions. 
Scattering polarised leptons off polarised nucleons pro- 
vides information on the spin composition of the nucleon. 
The major aim of the HERMES experiment is to deter- 
mine the spin contributions of the various quark flavours 
to the spin of the nucleon from the combination of inclu- 
sive and semi-inclusive deep- inelastic polarised scattering 
data @]. 

At centre-of-mass energies where weak contributions 
can be neglected, inclusive polarised deep-inelastic scat- 
tering is characterised by two spin structure functions: 
gi{x, Q 2 ) and g2(x, Q 2 ). Here Q 2 is the negative squared 
four-momentum of the exchanged virtual photon with 
energy v and x = Q 2 /2Mv is the Bjorken scaling vari- 
able, where M is the nucleon mass. The fractional en- 
ergy transferred to the nucleon is given by y = v j ' E for a 
lepton beam energy E. In leading order QCD the struc- 
ture function <?i is given by the charge weighted sum over 
the polarised quark (anti-quark) spin distributions Aqj 
(A?/): 

9i(x,Q 2 ) = ^$>f (Ag^Q 2 ) + Aq f (x,Q 2 )). (1) 
/ 

Here e/ is the charge of the quark (anti-quark) of flavour 
/ in units of the elementary charge and x is interpreted as 
the fraction of the nucleon light-cone momentum carried 
by the struck quark. 

The focus of this letter is the proton structure function 
g\, which provides a powerful constraint on the polarised 
quark distributions. The spin structure functions g\ and 
<72 can be determined from measurements of cross section 
asymmetries by combining data taken with a longitudi- 
nally polarised lepton beam and different spin orienta- 
tions of the target nucleons. Scattering off a longitudi- 
nally polarised target or a transversely polarised target 
yields the asymmetries An and A±, respectively. This 
letter reports on the results for g^ from the measure- 
ment of ^4 1| using inclusive deep-inelastic scattering data 
collected in 1997. For this HERMES measurement the 
target was pure polarised hydrogen gas without dilution 



by other unpolarised atomic species, in contrast to pre- 
vious measurements of g\ performed with solid targets 
at SLAC § and at CERN The target thickness, 

fraction of polarisable material, and the polarimetry are 
all significantly different for solid and gaseous targets. 

The HERMES experiment is located in the East 
straight section of the HERA storage ring at the DESY 
laboratory in Hamburg. It uses the positron beam of 
27.57 GeV energy with beam currents decreasing typi- 
cally from 40 to 10 mA in eight hours. The positrons 
become transversely polarised by the emission of syn- 
chrotron radiation p). Longitudinal polarisation of the 
positron beam at the interaction point is achieved with 
spin rotators p] situated upstream and downstream of 
the HERMES experiment. Equilibrium polarisation val- 
ues in the range of 0.40 to 0.65 are reached with a rise- 
time of about 30 minutes. The beam polarisation is 
continuously measured using Compton back-scattering of 
circularly polarised laser light, achieving a statistical ac- 
curacy of typically 1% in 60 s. Two polarimeters are used, 
one measuring the transverse polarisation in the HERA 
West straight section |7j] and the other measuring the lon- 
gitudinal polarisation near the HERMES target M. Af- 
ter a normalisation by rise-time measurements they give 
consistent results for the entire running period The 
systematic uncertainties are respectively 3.4% and 4.3% 
of the measured values, dominated by the normalisation 
uncertainty of 3.3% as determined from the rise-time cal- 
ibration. Data were accepted for this g\ analysis when 
the polarisation value was above 0.30; the average polari- 
sation for this data was 0.55. For about 12% of the data, 
only the longitudinal polarimeter was operational. 

The HERMES polarised proton target is formed 

by injecting a nuclear-polarised beam of atomic hydro- 
gen from an atomic beam source (ABS) |L2| into a tubu- 
lar open-ended storage cell. The cell confines the atoms 
in the region of the circulating beam and increases the 
probability of a positron-proton interaction by a factor 
of approximately one hundred compared with the free 
atomic beam. The storage cell is made of 75 (im thick 
aluminium, 400 mm long, and has an elliptical cross sec- 
tion 9.8 mm high and 29 mm wide. It is cooled to 100 K 
giving an areal density of 7x 10 13 atoms/cm 2 . The proton 
polarisation in the atomic beam is above 0.95, while the 
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electron polarisation is less than 0.05. The beam also 
contains unpolarised hydrogen molecules at a level less 
than 1%. Due to depolarisation mechanisms which take 
place within the cell, the proton polarisation within the 
cell is slightly lower than in the atomic beam. These are 
spin exchange in atom-atom collisions, as well as depo- 
larisation and recombination to molecules during inter- 
actions with the cell walls. The cell walls are coated with 
drifilm to minimise the wall collision effects. There is 
good evidence that recombination is further suppressed 
by water adsorbed on the cell surface during normal oper- 
ation |H],[l4|] . A static magnetic field directed parallel to 
the positron beam axis is provided throughout the cell to 
define the quantisation axis. The operational field value 
is chosen to avoid resonant depolarisation of the protons 
by the pulsed magnetic field caused by the bunch struc- 
ture of the HERA beam jlB). The spin direction in the 
target can be reversed in less than one second by selecting 
different spin states in the ABS. In operation the length 
of the time interval between reversals was randomised 
and was of the order of one minute. The atomic beam 
is injected into the storage cell via a side tube connected 
to its centre. The atoms diffuse to the open ends of the 
cell generating a triangular density distribution along the 
beam axis. The escaping gas is removed from the storage 
ring by high speed vacuum pumps. The residual gas in 
this vacuum system produces a further small source of 
unpolarised hydrogen molecules in the storage cell. The 
gas at the centre of the cell is sampled via a second side 
tube. Both the nuclear and electron polarisations of the 
atoms in this sample are measured with a Breit-Rabi po- 
larimeter (BRP), and the atomic fraction with a target 
gas analyser (TGA) jll|] . These measurements are made 
continuously achieving a statistical accuracy of typically 
0.02 in 60 s. 

The target proton polarisation Pr as seen by the 
positron beam is given by 



P T = oi [a T + (1 - a T )0\ B 



■atom 
T 



(2) 



Here P| tom 



is the proton polarisation in the atoms; its 
value was 0.92 ± 0.03. The atomic fraction oto accounts 
for the presence of the small number of molecules origi- 
nating from the ABS and the residual gas in the vacuum 
system; it had a value of 0.99 ± 0.01. The fraction of 
atoms that are in the form of molecules produced by re- 
combination is (1 — a r ). The value of a r was 0.93 ± 0.04. 
The values of P|; tom and a r were corrected for the sam- 
pling efficiencies of the BRP and TGA. These corrections 
depend on the measured values and the knowledge of de- 
polarisation mechanisms inside the cell. The quantity (3 
is defined as the ratio of the polarisation of protons in 
molecules from recombination to the polarisation of pro- 
tons in the atoms. Limits on this ratio (0.2 < j3 < 1.0) 
were derived p4| from the relationship between the tar- 
get polarisation and the corresponding asymmetry Am 



measured when the recombination rate was high due to 
deliberately modified conditions at the cell surface. The 
resulting target polarisation Pr was 0.88 ± 0.04. The 
quoted values are averages over the data taking period 
and the uncertainties are systematic. 

The HERMES detector is a forward spectrometer with 
a dipole magnet providing an integrated field of 1.3 Tm. 
The magnet is divided into two identical sections by 
a horizontal iron plate that shields the positron and 
proton beams from the magnetic field. Consequently, 
the spectrometer consists of two identical detector sys- 
tems, and the minimum polar angle for the acceptance 
of scattered positrons is 40 mrad. The maximum angu- 
lar acceptances are ±140 mrad vertically and ±170 mrad 
horizontally. For tracking in each spectrometer half, 
42 drift chamber planes and 6 micro-strip gas chamber 
planes are used. Fast track reconstruction is achieved 
by a pattern-matching algorithm and a momentum look- 
up method |16|. For positrons with momenta between 



3.5 and 27 GeV, the average angular resolution is 0.6- 
0.3 mrad and the average momentum resolution Ap/p 
is 0.7-1.3% aside from bremsstrahlung tails. A detailed 
description of the spectrometer is found in Ref. |l7[] . 

The trigger is formed by a coincidence of signals from 
three hodoscope planes with those from a lead-glass 
calorimeter, requiring an energy of greater than 1.5 GeV 
to be deposited locally in the calorimeter. Positron iden- 
tification is accomplished using the calorimeter, a scin- 
tillator hodoscope preceded by two radiation lengths of 
lead, a six- module transition radiation detector, and a 
C4F10/N2 (70:30) gas threshold Cerenkov counter. 

An adjustable two-stage collimator system is mounted 
upstream of the target cell to protect target and spec- 
trometer from synchrotron radiation and from beam 
halo. The number of triggers originating from parti- 
cles scattered from the storage cell walls was negligi- 
ble. The luminosity is measured by detecting electron- 
positron pairs from Bhabha scattering off the target gas 
electrons, in two NaBi(W04)2 electromagnetic calorime- 
ters, which are mounted symmetrically on either side of 
the beam line. 

The cross section asymmetry A|| was determined using 
the formula 



An = 



N~L+ -N+L~ 



Lr 



N-L+ +AH 



(3) 



Here N + (N~) is the number of scattered positrons for 
target spin parallel (anti-parallel) to the beam spin orien- 
tation. The deadtime corrected luminosities for each tar- 
get spin state are and Lp, the latter being weighted 
by the product of the beam and target polarisation values 
for each spin state. The luminosities were corrected for 
a small Bhabha cross section asymmetry caused by the 
typically 3% polarisation of the electrons in the target. 
The effect on A\\ was determined to better than 0.2%. 
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The kinematic requirements imposed on the data were: 
Q 2 > 0.8 GeV 2 , 0.1 < y < 0.85, an energy of the hadronic 
final state W > 1.8 GeV, and a minimum calorimeter en- 
ergy deposition of 3.5 GeV. After applying data quality 
criteria, 1.7 x 10 6 events were available for the asymme- 
try analysis. The data were analysed in bins of x and 
also y in order to include the effect of the variation of the 
virtual photon depolarisation factor D with y. For each 
bin and spin state the number of scattered positrons was 
corrected for e + e~ background from charge symmetric 
processes. This correction was at most 8% in the bins 
with x < 0.08 and y > 0.6, and negligible for the remain- 
der of the kinematic range. For the average positron 
identification efficiency of 99%, the average contamina- 
tion of misidentified hadrons was negligible with values 
not greater than 1% in the lowest x bins. The asymme- 
try was further corrected for smearing effects due to the 
finite resolution of the spectrometer, which were deter- 
mined by Monte Carlo simulations to be about 8% at low 
x and in the range of 2-3% at high x. QED radiative cor- 
rection factors were calculated following the prescription 
given in Ref. The corrections were determined to 

be less than 8% at low x and high y, decreasing to less 
than 1% at higher x values. Both smearing and radiative 
corrections were calculated using an iterative procedure. 

The structure function ratio g\/F\ is approximately 
equal to the longitudinal virtual photon asymmetry A\. 
It was calculated in each (x, y) bin from the longitudi- 
nal asymmetry A\\, corrected for the effects mentioned 
above, using the relation 
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D 



+ (7 - n)M 



(4) 



The virtual photon depolarisation factor D = [1 — (1 — 
y)e]/(l+eR) depends on the ratio R = ox /or of longitudi- 
nal to transverse virtual photon absorption cross sections 
and is approximately equal to y. The kinematic factors 
are defined as e = [4(l-y)-7 2 y 2 ]/[2y 2 + 4(l-y)+7 2 y 2 ], 
7 = 2Mx/ y^Q 2 , and r\ = ejy/[l — e(l — y)}. The mag- 
nitude of the transverse virtual photon absorption asym- 
metry A\ has been measured previously to be small || . 
Its contribution to gf/Ff is further suppressed by the 
factor (7 — rj) and was taken into account using a fit 
based on existing data |2|,(§: A% = 0.5 • x/ \[Q 2 . 

The x dependence of the structure function ratio 
g\/Ff at the measured (Q 2 ) for each value of x is shown 
111 Fig. 0. The averaged values of x, Q 2 , and the struc- 
ture function ratio g\/Ff are listed in Table |[ Fig. || 
shows a comparison of the measured g\ / Ff values with 
recent results of previous experiments, E-143 0] and 
SMC Q]. There is good agreement between these three 
sets of data, although the (Q 2 ) values of E-143 and HER- 
MES differ from those of SMC by a factor between five 
and ten. This demonstrates that there is no statistically 
significant Q 2 dependence of the ratio g\/Ff in the Q 2 



range of these experiments. The spin structure function 
g^(x,Q 2 ) was extracted from the ratio g\/Ff using the 
relation F% = F^l + 7 2 )/(2x(l + R)) with parameterisa- 
tions of the unpolarised structure function F^{x, Q 2 ) |l9| ] 
and R(x,Q 2 ) [p0|. The values of g\ at the measured x 
and Q 2 and after evolution to a common Qq of 2.5 GeV 2 
are also given in Table ||. The evolution was done un- 
der the assumption that the ratio of g\/Ff does not de- 
pend on Q 2 . The same assumption was used to evolve 
the HERMES data to the Qq values as published by E- 
143 and SMC. As shown in Figs. ||a and ||b, the evolved 
HERMES data are in excellent agreement with both data 
sets. The apparent Q 2 dependence of g\ seen when com- 
paring Fig. |3^, with Fig. |^b entirely originates from the 
Q 2 dependence of the unpolarised structure functions Fj 5 
and R. 
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FIG. 1. The structure function ratio gf/Ff of the proton 
as a function of x, given for the measured (Q 2 ) at each value 
of x. The error bars show the statistical uncertainties and the 
band represents the total systematic uncertainties. 



The systematic uncertainty of this measurement of 
g\/Ff, illustrated by the band in Fig. 0, is about 8% 
over the entire x range. The dominant sources of system- 
atic uncertainties are the beam and target polarisation 
discussed above. An uncertainty of 2.5% is included to 
account for the uncertainty in the description of the spec- 
trometer geometry. The uncertainty originating from the 
combined effects of smearing and QED radiative correc- 
tions is between 1% and 4%. The data were searched for 
possible systematic fluctuations in the measured asym- 
metry. This was carried out by dividing the data into 
sub-samples defined by various parameters, e.g. time, 
detector geometry and beam current. Within the statis- 
tical accuracy of the data no additional systematic effect 
was detected. 



4 




FIG. 2. Comparison of the data shown in Fig. Qwith recent results for gf/Ff obtained at SLAC (E-143) and A? at CERN 
(SMC) for Q 2 > 1 GeV 2 . The inner error bars show the statistical uncertainties and the outer ones the quadratic sum of 
statistical and total systematic uncertainties. 
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FIG. 3. The spin structure function g\ of the proton as a function of x. The HERMES data are evolved to a) Q 2 , — 2 GeV 2 
and b) Q 2 , = 10 GeV 2 assuming gf/F^ to be independent of Q 2 . This measurement is compared to recent results for 
Q 2 > 1 GeV 2 from E-143 || and from SMC [§, the latter shown for x > 0.01 only. The error bars are defined as in 
Fig. | 



The effect of the assumption for on the extracted 
ratio gi/Ff was estimated from an analysis of the varia- 
tion within the statistical and systematic uncertainties of 
the recent A^ data of E-143 0], the assumption = 0, 
and the use of the ansatz for given in Ref. pl| . The 
systematic uncertainty due to this effect was estimated 
to be at most 1.5%. 



The values of the ratio g\jF\ and also g\ derived 
from Am depend on R. Using the uncertainties for R 
given in Ref. p^] changed the g\/F^ (g\) values by typ- 
ically 3% (2%) and at most 5.5% (3.6%) at x < 0.1. 
In the HERMES kinematic region this uncertainty range 
of R covers the R measurement recently published for 
x < 0.12 P^j , Variations in the i 7 ^ parameterisation 
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calculated in Ref. |19|] including the normalisation uncer- 
tainty of the F% values of 0.7% resulted in changes in the 
extracted g\ of about 2.5%. The influence of R and 



has not been included in the systematic uncertainty of 
the measurement and the values are given separately for 
each x bin in Table Q. 



TABLE I. Results on gf(x,Q 2 )/Fl(x,Q 2 ), gf(x,Q 2 ) and gl{x,Q 2 ) evolved to Ql = 2.5 GeV 2 assuming gf/Ff to be 
independent of Q 2 . The measured (Q 2 ) values are given in GeV 2 . For each x bin the uncertainties concerning R and F2, S(R) 
or 8(R,F2), are listed separately. The systematic uncertainties are dominated by those of the beam and target polarisation 
which constitute a normalisation uncertainty only. 



(x) (Q 2 ) g{ /F 1 p zbstat.zbsyst.zb(5(.R) g P zbstat.zbsyst.zb<?(i?, F 2 ) g\{x, Qjj)±stat.±syst.±5(fl, F 2 ) 
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Nucleon spin structure functions may be characterised 
in terms of sum rules which involve integrals of g\ over x 
for a given Qq p3j . The integral of g\ (x) evaluated in the 
measured region 0.021 < x < 0.85 and Q 2 > 0.8 GeV 2 
is 0.122 ± 0.003(stat.)±0.010(syst.) at Ql of 2.5 GeV 2 . 
In the integration the x dependence of Ff within the 
individual x bins was fully accounted for, whereas the x 
dependence of g\ / Ff was treated in first order only. Pos- 
sible Q 2 dependences not excluded by the present statis- 
tical accuracy were investigated using Q 2 dependent fits 
to the gi/Ff data according to Ref. (|] and a Q 2 de- 
pendence of as parameterised by a next-to-leading 
order QCD analysis [Q. The effect on g\ was at most 
8% at the lowest x and Q 2 values and was included in 
the systematic uncertainty of the integral. A contribu- 
tion to the systematic uncertainty of 0.003 arose from 
the uncertainties due to R and described above. The 
result obtained for the integral was compared with those 
from E-143 @ and SMC f§, both calculated with the 
same integration scheme and for the kinematic range of 
HERMES. With respect to statistical uncertainties the 
agreement is better than 0.6a and 1.2<r, respectively, well 
inside the normalisation uncertainty given by each exper- 
iment. 



In summary, the proton structure function ratio g\ / Ff 
was measured with good statistical and systematic preci- 
sion. The results were obtained using an entirely differ- 
ent technique from that used in all previous experiments, 
involving a longitudinally polarised positron beam in a 
high energy storage ring and an internal polarised pure 
hydrogen gas target. The data are in good agreement 
with those obtained with solid targets at both similar 
and much higher Q 2 values, indicating that the system- 
atic uncertainties are well understood for both techniques 
over the entire measured x range. Evolved to the same 
values of Q 2 all recent data on g\ and on its integral are 
consistent. 
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